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AN APPROXIMATE SOLUTION FOR LUMINOSITY DISTANCE IN 


ZERO-PRESSURE RELATIVISTIC MODEL UNIVERSES 

THAT HAVE THE U-PROPERTY 

By Windsor L. Sherman 
Langley Research Center 


SUMMARY 


Finite-density zero-pressure models of the universe based on general rela- 
tivity are studied. These models have the U-property (time-dependent uni- 
formity) and a metric subspace described by the Robertson-Walker metric. Highly 
accurate approximate closed-form expressions for the radial metric variable, the 
luminosity distance, the time of light travel, and the scale factor of the "uni- 
verse are obtained. The results are used to obtain new forms of the redshift- 
magnitude and count -magnitude relations that are useful for the analysis of 
observational data once the signs of the cosmical constant and curvature con- 
stant are known. The redshift -magnitude relation is of a form that permits the 
determination of the density parameter from apparent magnitude and redshift data 


INTRODUCTION 


In order to find the model universe that best represents the observed uni- 
verse, it is necessary to determine certain parameters from observational data. 
These parameters include the acceleration and density parameters. At present 
the most satisfactory expression to use for the determination of these param- 
eters is the redshift-magnitude relation. Several forms of the redshift- 
magnitude relation have been derived and are presented in references 1 to 3* 

Two of these redshift-magnitude relations are based, either explicitly as in 
reference 1 or implicitly as in reference 2, on the assumption that a term 
called the cosmical constant is zero. Reference 3 presents a redshift-magnitude 
relation based on a model universe that has a nonzero cosmical constant. How- 
ever, because of the assumption of zero material density used in deriving the 
model universes, some of the zero-density models presented in reference 3 are 
not suitable for the analysis of observational data. The unsuitable models 
occur, for the most part, when the acceleration parameter is less than or equal 
to zero. Because models with an acceleration parameter that is less than zero 
may be important in the analysis of observational data when evolutionary effects 
are considered, finite-density zero-pressure model universes that possess time- 
dependent uniformity (U-property) are studied to determine whether the catalog 
of redshift-magnitude relations for model universes with a nonzero cosmical 
constant could be improved. 



The study resulted in a very good approximate expression for luminosity 
distance, and hence, for the redshift-magnitude relation, for model universes 
with a nonzero cosmical constant and finite density. The only model universes 
not included are those with a zero cosmical constant. For the model universe 
with a zero cosmical constant an exact solution exists. (See ref. 1.) 


ANALYSIS 


Distance in the observed uni verse is given by 


m - K = 5 log 1Q D ? + £M 0 + M 0 - 5 


( 1 ) 


where m is the apparent magnitude of the source, K is the correction for 
redshift to apparent magnitude, M Q is the absolute magnitude of the equivalent 

local source, Z^Mq is the evolutionary correction to the apparent absolute mag- 
nitude, and Dj is the luminosity distance. (See appendix for a complete sym- 
bol list.) The luminosity distance is obtained from the model universe, which, 
in the present study, is based on general relativity. In relativistic model 
universes (ref. 2), is given by 

D z = Rq(1 + 5)s(cd) (2) 

where R c is the present value of the scale factor, 5 is the redshift, co is 
a function of the radial metric variable (ref. 4), and S(cd) is a function of 
co that depends on the curvature of space. If the model also possesses zero 
pressure, time-dependent uniformity, and a metric subspace defined by the 
Robertson-Walker metric, S(co) is given by 

S(cd) = -i- sinh /-k a) (3) 

fk 

where k is the curvature constant. In this case the variable co is given by 
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where c is the speed of light in a vacuum, cr Q is the density parameter, H 0 
is the Hubble parameter, and A is the cosmical constant. The constants 

Ro 2 

and A are given by 
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A = 3H d 2 (o q - Qq) (6) 

where q Q is the acceleration parameter. Equations ( 5 ) and (6) are used to 
v 2 

eliminate and A from the integrand of equation (4), and the following 


equation results: 
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(7) 


This integral is elliptic and cannot he integrated in terms of simple func- 
tions. Mattig (ref. l) simplified this integral by assuming that A = 0 and 
thus obtained simple closed- form solutions. In reference 3 the integral was 
simplified by assuming that the density parameter a Q was zero and showed that 
under certain conditions the resulting zero-density models could be applied to 
the observable universe. The assumption that A = 0 deletes the constant 
term in equation ( 7 ), whereas the zero-density assumption affects terms of all 
orders, including the constant term. As shown in reference 3* there are cer- 
tain zero-density models, primarily those with 0 ^ % -1, that cannot be 

used for the analysis of observational data. Because of this difficulty, equa- 
tion ( 7 ) was studied to determine whether simple closed-form solutions that were 
good approximations of equation (7) could be obtained without setting the den- 
sity equal to zero. 


There are several approaches for obtaining approximate solutions in terms 
of simple functions. One of the approaches is to expand the integrand in a 
MacLaurin Series about 6=0, and another is to determine whether higher order 
terms involving density, such as 2a y can be neglected. The second approach 
is used in the present analysis. 

If the term 2a 0 (l + 5)5 in equation ( 7 ) could be neglected, S(od), and 
consequently the luminosity distance, would be algebraic expressions. An anal- 
ysis of the relative significance of the terms under the radical in equation (7)> 
however, showed that the term 2a 0 (l +5)5 was significant, and this approach 
was abandoned. 


The denominator of equation ( 7 ) was expanded so that 5 became the varia- 
ble. The expansion gave 


CJD = 


An analysis of the relative significance of the terms in the denominator of 
equation (7a) showed that for the presently accepted ranges of a Q and q^ 


H 0 R 0 


/ 

v-/ r. 


d6 


[2a 0 6 5 + (3 a Q + + l) 5 2 + 2 (q Q + l) 6 + l] 


1/2 


(7a) 


3 



(0.016 ^ o Q ^ 0.16 and 0.5 = % ^ 2. 5); the term 2a 0 s5 made a negligible 
contribution to the denominator when 8 < 1. 


The general solution of equation (7a) with 2a 0 &3 =0 is 

-1 (1 + q Q + 3 ct 0 )5 + (q Q + 1) 
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where the coefficient of the integral has been written as 


£ l - 3 g o 
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through the use of equation (5). In figure 1 is a comparison of — given 



by equations (7a) and (8) for the 
model universe for which A < 0 
and k = -1. For both models shown 
a 0 = 0.08; for one model q Q = 0.5* 
whereas for the other , q 0 = 2.5. 

As shown in figure 1 , differences 
between the exact and approximate 
forms cannot be detected for 
q 0 = 2.5j and they are barely 
detectable for q 0 = 0.5 at the 
higher values of 8. This fact is 
interesting because the term 2a 0 83 
was neglected on the basis that 
8 < 1. The results shown in fig- 
ure 1 indicate that, for present 
conditions, the term 2a 0 &3 is 
negligible out to 8=1, and the 
results of calculations (not shown) 
indicate that this condition 
extends appreciably beyond 8=1. 
This occurrence greatly extends the 
usefulness of the solution for to. 

The function S(u) is 
obtained by substituting equa- 
tion (8) into equation (5) and is 
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Figure L- Comparison of integration of equation 7(a) with integration 
of equation (8) for c 0 = 0.08. 
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Because, as shown in figure 1, the results for cd in the exact and approximate 
solutions are very nearly equal, there should be little difference in S(cjd) for 
the exact and approximate expressions because the hyperbolic sine is operating 
on arguments that are almost equal for the exact and approximate solutions . 


The luminosity distance is obtained by substituting equation ( 9 ) into equa- 
tion (2) and is as follows: 
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In figure 2 is compared the approximate solution for D 1 which is given by equa- 
tion (10) with the exact solution for which is obtained by using the value 

of cd given by equation (7a) in equations (2) and ( 3 ). As can be seen, the 

cdH 0 R 0 

agreement is as good as that found for — 

This approximation is good out to 5 = 1, as 
shown by figures 1 and 2, in spite of the fact 
that it was obtained under the restriction 
that 5 < 1. Unfortunately, this approximate 
solution holds only when A ^ Oj however, 
Mattig’s work (ref. l) provides an exact solu- 
tion for A = 0. 



Figure 2.- Comparison of exact and approximate luminosity 
distances. o 0 = 0.08; H 0 = 100 km/sec/Mpc. 


Even though the solution for given 

by equation (10) is a very good approximation 
to the exact solution for D^, differences are 
introduced as 5 increases above 1. Let € 
be the difference between the exact and 
approximate forms of equation ( 9 ); then D^ e , 
the difference in 
e, is given by 
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where <x> is given by equation (8). The difference e could be read from a 
plot similar to figure 2 extended to higher values of Sj thus at the higher 
values of S the deleted cubic term might have introduced a difference that 
does not show in figure 2. In reference 3 the difference between the values of 
R 

— 2-2- for the finite- and zero-density models was 0.005 at 5=1. This dif- 
ference, which corresponded to a difference in D j of 10® parsecs, was con- 
sidered acceptable. For equation (ll) D^ € was assumed to be 1C)8 parsecs, 
and € was assumed to be 0.005- For these values and for a Q = 0.08, equa- 
tion (ll) indicated that the errors would occur at 8 = 2.25 and 8 = 4.95 
for q 0 =0.5 and q 0 = 2.5j respectively. In these calculations H 0 was 
taken as lOO/km/sec/Mpc . 

Equation (10) represents a good closed-form approximate solution for the 
luminosity distance in a zero-pressure finite-density model universe that pos- 
sesses the U-property and has a nonzero cosmical constant. Unfortunately, there 
are many model universes covered by this equation. (See ref. 5*) Each of these 
universes involves a slightly different form of equation (10). In addition, 
singularities arise which must be handled as indeterminate forms. If the signs 
of k and A are known, equation (10) is useful in establishing the model uni- 
verse, because only one specific form of equation (lO) is necessary. 

In addition to the luminosity distance, the time of light travel and the 
scale factor R are of interest . The fact that the term 2a Q 83 can be neg- 
lected makes it possible to derive accurate closed-form expressions for the 
variables t Q - t and R. 

The integral for t Q - t from reference 3 is 
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where t 0 is the present time and t is some other time in the past. When the 
term 2a 0 83 is neglected, the expression integrates to 
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which is the general solution for t Q - t in the model universe with a nonzero 
cosmical constant . 


Figure 3 compares the numerical solution of the integral (eq. (l2)) with 
the approximate solution given by equation (l2a) for k = -1 and A < 0. In 
this figure (t Q - t) H q is plotted as a function of 8, and, as can be seen, 
the agreement between equations (l2) and (l2a) is excellent. 
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Figure 3.- Comparison of solutions for t 0 - ^Hq" 1 given by elliptic 
integral equation (11) and the approximate solution, equation (12). 


The inversion of equation (12a) gives a solution for R, which is 
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RELATIONSHIPS CONNECTING- OBSERVATION AND THEORY 


(13) 


The two most useful relationships connecting observation and theory are the 
redshift-magnitude relation and the count -magnitude relation. Of the two, the 
redshift -magnitude relation is more useful with present data and provides the 
stronger connection between observation and theory. 

The redshift-magnitude relation is obtained by substituting equation (10) 
into equation (l): 
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M 

To the first approximation is given by - 6, M Q being the change in 

absolute magnitude per epoch due to evolutionary effects. (See ref. 2.) This 
form of the redshift -magnitude relation contains three 'undetermined constants, 
q 0 , a 0 , and H Q , and can, in theory, be used with observational data in a 

least-squares computing process to find the best values of these parameters. 

Figure 4 shows the effect of varying q 0 when Hq = 100 km/sec/Mpc, 
o Q = 0.08, and M Q = 20. 56 . The separation of these curves varies from 0.127 m 
between the curves for q Q = 0.5 and q Q = 1.0 to 0.078 m between the curves 
for q Q = 2.5 and q 0 = 3*0 at 6 = 1. This separation is slightly smaller 
than the separation for the zero-density model universe (ref. 3)> and indicates 
that either observation must be extended beyond 6=1 or the precision must 
be improved before a model of the universe can be defined. The separations are 
still less than the best root mean square of the residuals found in reference 3* 
The calculations in reference 3 are considered to be applicable to this case 
because there is little or no effect of a 0 on the separation curve out to 
redshift values of 0.46. 



m- K 

Figured- Effect of varying q 0 on redshift-magnitude relation (eq. (14)). H 0 = 100 km /sec/Mpe; o 0 = 0.08; M 0 = 20.56; AM = 0, 
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Figure 5 shows the effect of cr Q on the redshift -magnitude relation when 
q Q = 1.5. The curves for cr Q = 0.016 and cr 0 = 0.l6 approximately correspond 
to the limits of Oort’s density range (see ref. 5) and the separation is about 
0.07 m at 5 = 1. Increasing a Q to 0.80 gave a separation of 0.286 111 at 
5=1. These results indicate that, up to a Q = 0.016, density has little 
effect on the redshift -magnitude relation. This conclusion confirms the results 
of reference 3 which indicate that up to o Q = 0.l6, the finite- and zero- 
density redshift -magnitude relations are the same. The curve for a Q = 0.l6 is 
separated very little from the curves for cr 0 - 0.0^ and o Q = 0.016 in fig- 
ure 5- This separation is not evident in reference 3 because of the scale used 
in plotting the figure, and for practical purposes it is negligible. It thus 
appears that the value of o Q must exceed 0.l6 before the effect of density 
becomes appreciable. 



m-K 

Figure 5.- Effect of varying o 0 on redshift- magnitude relation (eq. (14)). H 0 = 100 km/sec/Mpc; q 0 = 1.5; M 0 = -20.56; AM = 0. 
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Another relationship that is useful in the analysis of observational data 
and that comes from the model universe is the count -magnitude relation. In the 
past, the count-magnitude relation has not been very satisfactory for the analy- 
sis of observational data because the distance surveyed by the 200-inch tele- 
scope is not large enough to give an adequate data sample. (See ref. 6.) How- 
ever, with the advent of the space telescope, this situation may change. 

The count -magnitude relation (ref. 3) is 


W(m) 


2jmR 0 ^ 

(-k)Qo 


^sinh fka) cosh 


CD 



(15) 


where CD is given by equation (8). In equation (l5), N(m) is the number of 
galaxies brighter than a given apparent magnitude m, n is the number of gal- 
axies per unit volume, and Q 0 is the number of square degrees in the celes- 
tial sphere. 


This count -magnitude relation for the finite -density model, unlike the one 
for the zero-density model (ref. 3), is very close to the exact count -magnitude 
relation because of the agreement of the exact and approximate values of cd 
and R 0 . 


CONCLUDING REMARKS 


Relativistic finite-density zero-pressure models of the universe that 
possess the U-property and have a metric subspace that is described by the 
Robert son- Walker metric have been studied. Approximate expressions for the 
radial metric variable and the luminosity distance have been derived. Up to 
redshifts of 1, the approximate solutions do not differ from the exact solu- 
tion obtained through the use of a high-speed digital computer. An expression 
for the difference in the luminosity distances that is valid for any redshift 
was obtained. In addition, similarly precise closed-form solutions have been 
obtained for the travel time of light from galaxies and the scale factor of 
the universe . 

These results are used to obtain a redshift -magnitude relation that holds 
for all finite-density zero-pressure model universes that have the U-property 
except for those with a zero cosmical constant. This expression is difficult 
to use because of the many models the solution contains and because of the sin- 
gularities that arise . However, if the signs of the cosmical constant and the 
curvature constant are known, the expression becomes very useful for the analy- 
sis of observational data. This form of the redshift-magnitude relation permits 
the calculation of the density parameter from apparent magnitude and redshift 
data. Lastly, a potentially useful form of the count -magnitude relation is also 
derived. 


Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., June 23, 19&5* 
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APPENDIX 

SYMBOLS 


Ho 

K 

k 

Mo 

*0 


m 

N(m) 

n 

% 

R 

Ho 

S(cd) 

t 

t O 

5 

A 

*o 

CD 


speed of light in vacuum 
luminosity distance 
Hubble parameter 

redshift correction to apparent magnitude 
curvature constant 
absolute magnitude 

evolutionary correction to absolute magnitude 


rate of change of absolute magnitude for present epoch, 

\ dt / t=t 0 

apparent magnitude 

number of galaxies brighter than apparent magnitude m 

number of galaxies per unit volume 

number of square degrees in celestial sphere 

acceleration parameter 

scale factor 

present value of scale factor 

function of o> that depends on curvature of space 
time before present 
present time 
redshift 

cosmical constant 
density parameter 

function of radial metric variable 
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